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ABSTRACT 

Here we present a number of improvements to weak lensing 3D power spectrum anal- 
ysis, 3D cosmic shear, that uses the shape and redshift information of every galaxy 
to constrain cosmological parameters. We show how photometric redshift probability 
distributions for individual galaxies can be directly included in this statistic with no 
averaging. We also include the Limber approximation, considerably simplifying full 3D 
cosmic shear analysis, and we investigate its range of applicability. Finally we show 
the relationship between weak lensing tomography and the 3D cosmic shear field itself; 
the steps connecting them being the Limber approximation, a harmonic-space trans- 
form and a discretisation in wavenumber. Each method has its advantages: 3D cosmic 
shear analysis allows straightforward inclusion of all relevant modes, thus ensuring 
minimum error bars, and direct control of the range of physical wavenumbers probed, 
to avoid the uncertain highly nonlinear regime. On the other hand, tomography is 
more convenient for checking systematics through direct investigation of the redshift 
dependence of the signal. Finally, for tomography, we suggest that the angular modes 
probed should be redshift-dependent, to recover some of the 3D advantages. 
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1 INTRODUCTION 



Cosmology is faced with a standard model that contains 
two dominant unknown components (dark energy and dark 
matter) and two untested assumptions (that general relativ- 
ity applies on large scales and that there was an inflation- 
ary period in the early Universe). In order to move forward 
from this unsatisfactory situation the need for high preci- 
sion methods, which exhibit a strong statistical signal to 
cosmological parameters (those that describe the standard 
model or deviations from it) to allow various models to be 
compared with a high discriminatory power, is obvious. 

It has become clear that 3D cosmic shear (Heavens, 2003; 
Castro, Heavens, Hitching, 2005; Heavens et al., 2006; Hitch- 
ing et al., 2007) and weak lensing tomography (Hu 1999; 
Amara & Refregier, 2007; Bernstein & Jain, 2006), in which 
galaxy redshifts and the weak lensing shear distortion are 
used simultaneously in the signal, is an approach which is re- 
markably sensitive to cosmological parameters through both 
the growth of structure and the geometric behaviour of the 
lensing effect itself. For example neutrino mass (Hannestaad 
et al., 2007; Hitching, Heavens, Verde, 2007; de Bernardis et 
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al., 2009; Jimenez et al., 2010), modified gravity (Thomas 
et al., 2008; Heavens, Hitching, Verde, 2008), dark matter 
(Camera et al., 2010; see Massey et al., 2010 for a recent 
review). In particular, given a current focus on dark energy 
(Albrecht et al., 2006; Peacock & Schneider, 2006) as the 
point in which we may find a departure from our standard 
model, it has been shown (and is widely accepted) that 3D 
cosmic shear has the potential to place tight constraints on 
the equation of state of dark energy. 

This article is concerned with the techniques outlined in 
Heavens (2003), Heavens et al. (2006) and Hitching et al. 
(2007) whereby the 3D shear field (a function of redshift 
and angle) is treated in its entirety. In this approach the 3D 
shear power spectrum is reconstructed using spherical har- 
monics and the cosmological signal is found in the covariance 
of the coefficients. This is an alternative to weak lensing to- 
mography where the 3D shear field is binned and projected 
in redshift space to produce a series of 2D projections, from 
which the auto- and cross-power spectra can be calculated 
(e.g. Hu 1999; Amara & Refregier, 2007, Schrabback et al., 
2009). 

3D cosmic shear has some advantage over tomography in 
that no binning is required, avoiding its consequent loss of 
information. One of the purposes of this article is to explore 
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the relationship between 3D cosmic shear and tomography. 
3D cosmic shear in principle allows a full sky reconstruction 
including the effects of sky curvature. Individual pairs of 
galaxies in redshift can be excluded from the estimator di- 
rectly (to reduce intrinsic alignment contamination; King & 
Schneider, 20021; Heymans & Heavens, 2003) and since ev- 
ery galaxy contributes directly and individually to the signal 
information such as the probability distribution in redshift 
of individual galaxies can be included in the estimator. 

This final point, the inclusion of individual redshift proba- 
bilities, has been alluded to but has not been implemented 
up until now in 3D cosmic shear. Previously the inclusion 
of photometric redshift uncertainties has been through the 
use of the average photometric redshift error distribution at 
each redshift. We will show that this approach can lead to 
a mis-estimate of cosmological errors, and introduces small 
biases in cosmological parameter. 

We will also present a Limber approximation of the full 3D 
shear field and show that for the range of (radial and az- 
imuthal) modes which are sensitive to cosmological infor- 
mation the approximation is very good. This approximation 
vastly simplifies the theoretical computation of the 3D co- 
variance, making it as simple as tomography to apply. We 
also explicitly derive the weak lensing tomography power 
spectra from the full 3D shear field, showing that these are 
related by two approximations (the Limber approximation, 
and the binning approximation), and demonstrate that to- 
mography provides a sampling of the physical wavenumber 
space which is £- and bin-dependent. 

This paper is organised as follows. In Section 2 we will review 
the 3D cosmic shear formalism and show how photometric 
redshift information from individual galaxies can be explic- 
itly included, using Fisher matrices we will investigate the 
impact on parameter error and bias. In Sections 3 and 4 we 
will present several approximations of the 3D shear field, us- 
ing the Limber approximation. We discuss results in Section 
5. 



2 PHOTOMETRIC 3D SHEAR ESTIMATOR 

3D cosmic shear decomposes the 3D shear field (a function 
of angular and radial position) using spherical harmonics 
(see Heavens, 2003; Castro, Heavens, Kitching, 2006; Heav- 
ens et al., 2006). Using the flat sky approximation we replace 
spherical harmonics with exponential functions, and assum- 
ing the Universe is flat we replace the more general ultra- 
spherical Bessel functions by spherical Bessel functions, a 3D 
shear estimator can be written, by summing over galaxies 



7i(M) = 



\ZIe* 



je(kr 



w{4) 



(i) 



from the photometric redshift assuming a fiducial cosmology. 
W(r) is an optional weight which we set to W — 1 for clar- 
ity in the remainder of this article, we note here that there 
may exist optimal weights, particularly for constraining a 
certain cosmological parameter set. r^ is a distance, not a 
redshift, so to convert from some data (spectra/photometry 
— > redshifts) requires the assumption of a cosmology - how- 
ever as shown in Kitching et al. (2007) this assumption is 
benign. Note that the transform convention of Castro et al. 
(2005) includes an extra factor of k, which we omit. The 
result of this expansion is a set of four 3D data vectors that 
are functions of angular wavenumber (£ — 2n/0) and ra- 
dial wavenumber and a real and imaginary part of 71 and 
72. Note that the behaviour of the Bessel functions (that 
je(x) <g; 1 for x <C €) ensures that a mode of radial wavenum- 
ber k also has no transverse power with smaller wavelengths, 
so we are justified in identifying k as a physical wavenum- 
ber in 3D space. The assumption of a flat sky (resulting in 
an exponential in place of 17™ above) can be relaxed (see 
Castro et al., 2005), resulting in a covariance of this esti- 
mator that is similar to the one that we derive below, wc 
leave a full investigation of this for a future article. The as- 
sumption of a flat Universe (resulting in Bessel functions in 
place of ultra-spherical Bessel functions) can also be relaxed, 
for computational and notational convenience we use Bessel 
functions in this article and take small perturbations about 
this model. 

Since the mean of the estimators is zero we use the covari- 
ance (Tegmark, Taylor & Heavens, 1997) of the estimator to 
constrain cosmology. In Appendix A we derive an expression 
for the expected covariance of the 3D shear estimators keep- 
ing track of any sum over galaxies. We may write the com- 
ponents of the shear estimator by 71 (k,£) — i (fn — I22) 

and 72 (fc,^) = Ti2 = T21, where the F terms are defined in 
terms of the lensing potential <f> by F^ = didj(j>, and have 
covariance 

c?/J(*i,fe,) = (t ij (k 1 ,e)t* j (k 2j £)) = 
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where 7? is the i th shear component for each galaxy g at 
angular position 9 and radial position r 9 ,, which is obtained 



F k (r,r') = S k (r-r')/[S k (r)S k (r')] (S k = sinh(fc), fc,sin(fc) 
for open, flat or closed geometries so F k (r, r') — (r — r')/rr' 
for a flat universe), a(r) — 1/[1 + z(r)] is the scale factor 
and P(k; r) is the 3D matter power spectrum. The prefac- 
tors to this expression are AQ the dimensionless area of the 
survey, A — 3fi m i/ 2 /2 where Q. m is the present dimension- 
less matter density and H is the present Hubble expansion 
rate. The £ 2 £ 2 is a factor that is introduced as a result of 
the derivative taken when transforming the lensing potential 
to the shear (see Heavens et al., 2006 and Kitching et al., 
2007). The I modes are uncorrelated due to the assumption 
of isotropy. To fully describe a non-Gaussian random field 
higher order potential correlations are required (see Munshi 
et al., 2010). The Tq and rj are calculated from photometric 
redshifts assuming a fiducial model, where the r 9 and r h are 
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Figure 1. Left panel: a sample of 50 redshift distributions from the sample used (from Bordoloi et al., 2009). Right panel: the n(z p ) = 
y"] p(z p ) from the total set of p(z), the dashed line shows the functional form for n(z) oc z 2 exp[— (1.4z/z m ) 1,5 ] from Taylor et al., (2007) 
with z m = 1.0 for comparison. 



their true values. This expression is one of the main results 
of this article. 

In making any parameter error forecasts we will use the 
Fisher matrix formalism (Tegmark, Taylor, Heavens, 1997) 
and assume Gaussian likelihood surfaces, and Gaussian dis- 
tributions for the data. For 3D cosmic shear the Fisher ma- 
trix was introduced in Heavens (2003) and examined in more 
detail in Heavens et al., (2006) and Kitching et al., (2007). 
Since we are assuming that the parameters affect the covari- 
ance, not the mean, the Fisher matrix is given by 



tp 9 



d<f> e J demic^Ct^Ce 1 



Ce,p] 



(3) 



where we include an integral over ^-space which includes a 
density of states in ^-space, g =Area/(27r) 2 (see Appendix 
B of Kitching et al., 2007). A comma represents a deriva- 
tive with respect to parameter a or /?. The covariance C 
is a contribution of signal terms and (shot) noise terms 

(J = (^signal + ^noise^ where ^ gigna , j g defined m equat ion 

(2) (we have suppressed the notation Cffy(ki,k2) — > Ci for 
clarity in equation 3) . The shot noise covariance matrix used 
is the same as that in Heavens et al., (2006). 



Throughout we will use the parameter set (with fiducial 
values) : fi m (0.3), fi de (0.7), iu (-0.95), w a (0.0), h(0.71), 
Sl 6 (0.045), <rg(0.8) and n s (1.0). We use the expansion of the 
dark energy equation of state as introduced in Chevallier 
& Polarski (2001) and Linder (2003), we include generally 



The density of states accounts for correlations between modes 
arising from partial sky coverage, equivalent to the / s fc H approach 
of many papers. Note that the insensitivity to large-scale modes, 
which is also a consequence of using a patch of sky, needs to be 
treated by a cut on £. The Fisher matrix approach assumes the 
data are Gaussian; see Munshi et al., (2010) for an investigation 
of non-Gaussianity in 3D shear. 



curved geometries where fit = I — fl m — fide 2 ■ Our matter 
power spectrum is constructed using the Eisenstein & Hu 
(1999) transfer functions, with baryon wiggles, and employ 
the Smith et al. (2003) non-linear correction; for curved cos- 
mologies with varying w(z) we use the iCosmo (Refregier et 
al., 2008) interpolation scheme (also discussed in Schrabback 
et al., 2009). 

Throughout we do not include any systematics effects in 
the predicted parameters errors. However we note that the 
impact from all primary systematics (photometric redshift 
calibration, intrinsic alignments and shape measurement er- 
ror) should be at most a factor of y/2 on cosmological pa- 
rameters (Kitching et al., 2008) in a realistic self-calibration 
regime using parameterised systematic functions. Using non- 
parametric approaches may degrade results further, but as 
shown in Kitching & Taylor (2010) and Kitching et al. 
(2009), using tomographic analyses, the overall cosmolog- 
ical impact is not severe. 

We assume a fiducial survey configuration of 20,000 square 
degrees, with a surface number density of 35 galaxies per 
square arcminute and an intrinsic ellipticity dispersion of 
0.3. Where we do not use direct photometric redshift proba- 
bilities we will use a redshift distribution of galaxies n(z) oc 
z 2 ex.p[-(lAz/z m ) 1 - 5 } given in Taylor et al., (2007) with a 
median redshift of z m = 1 and a Gaussian redshift disper- 
sion with a redshift error a z (z) = 0.03(1 + z). This survey 
configuration is similar to the proposed Euclid wide survey 
(Refregier et al., 2010). Throughout we use a maximum az- 
imuthal wavenumber of ^ max = 5000 and a maximum radial 
wavenumber of fc ma x = 1.5hMpc _1 to avoid the highly non- 
linear regime where theoretical predictions for the power 
spectrum may be unsound (e.g. Rudd et al., 2008; Joudaki 
et al., 2009). 



2 In this case ultra spherical Bessel functions should be used, 
but as pointed out in Kitching et al. 2007 in the t 2> 1 and 
k 2> (curvature) -1 regime the use of the ji(kr) is well justified. 
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2.1 Including Photometric Redshifts 

For the depth of surveys required for optical weak lensing, 
it is impractical to obtain spectroscopic redshifts, so photo- 
metric redshifts z p are obtained, typically from broad-band 
photometry. Some photo-z codes return a complete posterior 
probability distribution for the redshift, given the photom- 
etry, and the purpose of this section is to show that this 
individual information can be used when it exists. Since it 
uses all the information available, it should yield a more 
accurate covariance matrix than the alternative of reduc- 
ing the redshift distributions to a simpler p(z\z p ); the form 
of p may be arbitrary (and include for example outliers), 
but its distribution is assumed in the alternative approach 
to be the same for all galaxies with the same photometric 
redshift. Note that, although a full posterior redshift distri- 
bution may be returned, we still do have to identify a unique 
photo-z for each galaxy, to provide an r% to insert into the 
spherical Bessel function weighting. 

To include galaxy photometric redshift errors we integrate 
over the posterior redshift distribution for each galaxy. We 
write this as p g (z\z g ) as a shorthand for the redshift prob- 
ability distribution given all the photometric information 
available for galaxy g: 

C?/J(fci,fei) = AMJ^-^A^Y. [^(*i' - o)J/(*2'-o)] 

g h 

dz'p g (z'\zg) / dz"ph(z"\z h ) 



df I df Fk{t ,f )Fk{t ,r ) 
Jo 

^ a(r'Mn Mk ' nMk ' f '' WP{k '' f ' )P{k "' f ' % (4) 



U e (r{z'),k') 



' "'df ^yw op^f) (5) 

a{r) 



where A 2 = AQ£ 2 £ 2 £ -t, A 2 . Equation (5) is a discrete case 
of the equations in Section 2.5 of Heavens et al. (2006). The 
ij refer to 7 combinations (see Appendix A), not redshift 
bins - this is a continuous 3D estimator. 



If we do not have individual posterior distributions, or wish 
to ignore the individual galaxy redshift error distributions, 
we can simplify p g (z\z g ) to a global p(z\z p ) so all galaxies 
at fixed photometric redshift are assumed to have the same 
distribution of true redshifts. In this regime the the G matrix 
is modified only so that 

Gf(ki,k')— dz p dz'je(kir(zp))n{zp)p{z'\zp)Ue{r(z'),k') (6) 



where p(z\z p ) is the redshift probability distribution at pho- 
tometric redshift z p , and n(z p ) is the number density of 
galaxies as a function of redshift. Equation (6) is in agree- 
ment with Castro et al. (2003) and Heavens et al. (2006). 

By including the individual galaxies in equation (5) we auto- 
matically take into account all overlap between photometric 
redshift posteriors, all outliers in the sample and have the 
best estimate for the covariance of the data vector. Of course 
the outliers have only been properly accounted for only if the 
p(z\z p ) is a correct estimate of the probability of a galaxy 
being an outlier; we leave an investigation into the effect of 
errors on the p(z\z p ) for future study. 



where r' — r(z'). This is the covariance of the full 3D shear 
field constructed only from a sum over the individual tracer 
galaxy population, and including all uncertainty that may 
be present from photometric redshift estimates. 

In this formalism we can retain individual posteriors. Since 
we do not want the data vector to depend on the cosmologi- 
cal parameters (and consequently violating the conditions of 
the Fisher matrix analysis), we must assume a fiducial cos- 
mology to translate a photo-z to a distance rf y However as 
long as this choice, for each galaxy, is the same as those used 
in the actual sum over data (equation 1) then the estimator 
will be unbiased. 

For clarity we can re-express the covariance as series of ma- 
trices 

O*i,fc0 = A 2 j^G?(k 1 ,k')G?{k 2 ,k') 

Gf(ki,k') = y^ji(fcirg) / dz'p g (z'\zg)Ut(r(z'),k') 



dz' 



^2je(kir%)p g (z'\z g ) 



Ue(r(z'),k') 



2.2 Forecasted Impact 

To illustrate the effects, we use a sample of realistic pho- 
tometric redshift probabilities {p(z p )} from Bordoloi et al. 
(2009) , which arc a simulated set of redshift probability dis- 
tributions like those expected from Euclid. We show these 
in Figure 1. We use a representative sample of 3000 redshift 
distributions, for computational speed, which for a surface 
number density of no = 35 per square arcmin represents a 
sample of galaxies from around 100 square arcmin. This is 
low compared to the expected number of galaxies in future 
surveys, however we use this sample as representative for 
the investigation of the approximations we highlighted in 
Section 2.1 - in fact with low(er) number statistics we may 
expect a larger deviation than in reality because the error 
on the mean will be larger (a 1/iV effect). To match the 
total number to the surface density we multiply the signal 
covariance by an extra factor of (no/3000) 2 (and the noise 
by no/3000). 

We consider two cases 



(i) First case, where we use the full covariance including 
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posteriors for the true redshifts 3 , equations (5). 

(ii) Second case, where we make an approximation by 
summing the individual posteriors in photometric redshift 
intervals b, to find a global number density distribution 
as a function of redshift p(z\z p ) = ^ p g (z\z p ), equation 
(6). This is one particular method of smoothing the red- 
shift distribution, we would expect similar results with other 
smoothing kernels with a similar redshift scale. In this ex- 
ample we use 20 bins in z p between the redshifts 0-2 (a step 
in redshift of 0.1). 



We show the n(z) from the sample of photometric redshift 
posteriors we use in Figure 1. In Figure 2 we show the impact 
on the 3D cosmic shear power spectra. The power spectra 
are a function of £ and two physical wavenumbers, hence 
we show the plane (fci,/c2) for a series of £ modes, as well 
as the diagonal Cf D (k, k). It can be seen that using an ap- 
proximation of the redshift distribution can cause significant 
residuals in the power spectra, with residuals of order 0.01-1 
over all scales. These residuals are most prominent at low-£ 
and small radial scales. 

Using the Fisher matrix approach Figure 3 shows how the 
approximation of averaging the posteriors p(z\z p ) in pho- 
tometric redshift, can affect cosmological parameter errors. 
In general the errors are degraded by averaging the individ- 
ual galaxy posteriors p g (z\z g ) by a factor of 10-50%, and 
in some cases the parameter degeneracies are changed. This 
is a result of the individual p g (z\z g ) distributions explicitly 
including outlying and non-Gaussian behaviour on a galaxy- 
by-galaxy basis. 

In making the approximation of p(z\z p ) the best fit cos- 
mological parameter values may also be biased, as well as 
the errors being affected. In Appendix C we show how to 
calculate the bias in the Fisher matrix approximation for 
the case that the parameter dependency is in the covari- 
ance (this is a generalisation of the result of Knox et al., 
1998). We find that the approximation made in this exam- 
ple introduces a negligible bias on all cosmological param- 
eters : fi m (-3.0 x 10 -5 ), ft de (2.5 x 10~ 4 ), w (1.3 x 10~ 4 ), 
io„(-1.2xlO -4 ), h(5.2x 10~ 4 ), Q 6 (l.lx 10 -6 ), cr 8 (7.3x 10~ 5 ) 
and n s (2.5 x 10~ 4 ). We note that changes in the covariance 
only impact cosmological parameter biases if the change in- 
duced is similar to the effect of any cosmological parameter 
(see for example Kitching et al. 2009; form filling functions). 
In this case the changes in the covariance are large, but the 
changes are not similar to the effect of cosmological param- 
eter so the biases are small. The amplitude of the bias will 
depend on the exact distribution and form of the posterior 
redshift distributions, here we use an example to introduce 



3 In practice we have used redshift-binncd true photometric red- 
shift posteriors, using the resolution provided by Bordoloi et al. 
(2009), which was 5000 bins in redshift. This is sufficient to il- 
lustrate the method, since this is much larger than the number 
of p(z) availabcl. A binned description is not a limitation of this 
method - one could imagine a basis set decription with no explicit 
redshift binning. 



the methodology, we leave a more detailed study for a future 
article. 

For any given experiment, the exact degradation will be de- 
pendent on the exact galaxy survey and quality of photo- 
metric redshifts available. 

Note that in Heavens et al. (2006) we made a further approx- 
imation by assuming that the p(z\z p ) = n(z p )pGa,uss(z\z p ) , 
where n(z p ) is the number density of galaxies as a function 
of photometric redshift and PGa,uss(z\z p ) is the normalised 
photometric redshift distribution at each redshift which we 
assummed to be a Gaussian at each redshift with a variance 
<Tz{z p )- 



3 LIMBER APPROXIMATION OF THE 3D 
SHEAR FIELD 



We now investigate various approximations of the 3D shear 
field. From LoVerde & Afshordi (2008) the Limber approx- 
imation can be encapsulated in the following substitution 



lim je(kr) 



2(^) 



kr - 



£+■ 



(7) 



where 8 D is the Dirac delta function. This is substituted 
into any integral or sum that the Bessel function appears; 
for example equations (5) and (6) . We note that this does not 
explicitly convert a 3D wavevector into a 2D wavevector, k — 
(k x ,k y ,k z ) — > (k x ,k y ) (which is explicit in an alternative 
and complementary derivation of the Limber approximation; 
Kaiser, 1998), but that this transformation is implicit when 
an integral over the delta function is performed. 

By substituting this approximation into the matrices in 
equation (5) we can find an approximate expression of the 
covariance in the high I limit. In this case the U matrix is 
modified to 



Ue irnbcl {r(z'),k') 



\/2(€+i) o(r„) 



where r v = {£ + \)jk' < r(z'). 



In the individual redshift error case the G matrix is modified 
to 



Gf{k u k') 



H'+h) 



Y, dz'p g (z'\zV)UY™ hcl (r(z'),k') (9) 

{j:r»=[«+i]/tl} 

where each galaxy only contributes to a given combination 
of fci and £ set by its comoving distance r^ = [£ + \\/k\. 
We note that this equation in practice would require some 
binning in k and £ to obtain a finite number of galaxies in 
this sum. 
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Figure 2. The 3D cosmic shear power spectra for the case where we use individual galaxy redshift distributions (upper panel), and 
using an approximation of the photometric probability distribution (second panels). The colour scale is shown in the bar above. The 
third set of panels shows the diagonal C^(fc, k) of the full 3D power spectra for the full case (black lines) and for the approximation (red, 
dark grey lines). The bottom panels show the modulus of the difference in the diagonal elements. Rightwards of the cusp in the bottom 
panels the difference between the power spectra becomes negative. For a Euclid survey with Bordoloi et al., (2009) p(z). 
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Figure 3. The Fisher matrix 2-parameter 1-cr constraints for the case where individual galaxy p(z) are used and the approximation of 
using a summed, approximated global p(z\z p ). The central cross shows the fiducial values. This is for a Euclid survey (Refregier et al., 
2010) using p(z) from Bordoloi et al., (2009). 
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If a global error distribution is used, the G matrix is modified 
to 

GC, Limber/ 7 i /\ / \ 

, (ki,k ) — —n(z v ) 

47T 

(10) 

where r v — (£ + \)/ki and z v are the photometric redshifts. 
The extra factor of (4 s ) comes from the fact that the 
Bessel integral in equation (6) is over z not r. 

By substituting the Bessel function approximation directly 
into the shear estimator, equation (1), it can be seen that 
the Limber approximation leads to an equivalence between 
the 3D k and I values and the tomographic r and £ values. 
We will discuss the Limber approximation further in Section 
4. 



3.1 Convergence of the 3D Cosmic Shear 



In Figure 4 we show the effect of the Limber approximation 
on the 3D cosmic shear power spectrum. We find that the 
Limber approximation is a remarkably good approximation 
to the full calculation for scales £ > 100, with residuals of ~ 
10~ 4 over all radial and azimuthal scales. The break at k ~ 
^/r ma x in the power spectra at each £, caused by the Bessel 
function inequality je(kr < £) ~ 0, is reproduced through 
the inequality expressed after equation (8) . For larger scales 
I < 100, there is a larger effect on the power spectrum. 

In Figure 5 we show the effect of the Limber approximation 
on the expected cosmological errors. We find that in most 
parameter combinations some information is inevitably lost, 
through the largest scales being down weighted, the increase 
in errors is between 1-30% with an average increase of ~ 
10%. 

We conclude that the Limber approximation is adequate 
for forecasting purposes, but if computer time allows, it is 
preferable to use the full expressions in data analysis. 



4 TOMOGRAPHY FROM 3D COSMIC SHEAR 

We now approximate the 3D shear field further and show 
how in the discrete real-space limit the tomographic power 
spectra can be reproduced. There have been some implicit 
references to this derivation, for example in Hu (1999), here 
we will show explicitly how the 3D shear field is related to 
the tomographic power spectrum. 

Weak lensing tomography is a flavour of 3D weak lensing in 
which the angular and redshift information of each galaxy is 
used. The practical distinction between 3D cosmic shear and 
tomography is that tomography divides the redshift range 
into a series of bins and the 2D shear transform in each bin 



is constructed. The auto (in a single bin) and cross (between 
bins) power spectra are used to constrain cosmological pa- 
rameters. 

Given the expressions in equations (8) to (10) we can now 
derive the weak lensing tomographic power spectra directly 
from the 3D shear field. In Appendix B we show how the 
3D cosmic shear power spectrum using the Limber approx- 
imation can be written as 



C e (fcl, 



mlH 4 



dr 



fe») = 
P(e/r;r)W{n,r)W(r 2 ,r) 



4c 2 / a 2 (r) 

where we define a weight factor as 

W(ri,r) — [n(ri)ri]r / dr'p(r') 



(11) 



(12) 



where n = l/ki and p(r[z]) = p(r[2]|r[z p ]) is the redshift 
probability distribution, equivalent to p(z\z p ) in equation 
(6). We condense the notation here and in Appendix B to 
match the literature for the tomographic case. This is still a 
full 3D estimator where n = r(zi) and r-z = r(z?) can take 
any value, and in practice would be a sum over all galaxy 
pairs. This is a key result of this article, using only two inte- 
grals a full 3D shear power spectrum can be computed from 
the 3D matter power spectrum (as simple as the standard 
tomographic approximation). 

Inspection of the previous equations shows that they are the 
usual expressions for the (auto- or cross-) power spectrum 
of tomography, from which we see that, under the Limber 
approximation, tomography samples discrete sets of physical 
wavenumbers, in an ^-dependent way: for shells at distances 

n, k - £/n. 

In summary to convert from 3D cosmic shear to weak lensing 
tomography we see that the following three steps must be 
taken 



• The Limber approximation must be applied to the 
full 3D shear estimator. 

• Fourier transform the kernel integration from har- 
monic space to real space. 

• Discretisation of physical modes through k — £/ri. 

The second step is benign in that no information should be 
lost, however the first and third steps do result in informa- 
tion loss (see Section 3.1 and 4.1). Interestingly for a specific 
redshift bin (r) and a specific azimuthal €-mode the tomo- 
graphic approximation only probes a single physical fe-modc 
k — £/r from the full 3D shear field; in contrast in 3D cosmic 
shear we have control over the k and £ modes over the whole 
redshift range. 

Clearly by fixing the distances of the tomographic bin- 
ning, we lose some flexibility over the physical wavenumbers 
probed, so there is a risk that either not all useful modes 
are included (increasing statistical errors), or that, for the 
nearby shells, the physical wavenumber range sampled ex- 
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Figure 4. The 3D cosmic shear power spectra for the full 3D cosmic shear calculation (upper panels), and using the Limber approximation 
(second panels). The colour scale is shown in the bar above. The third set of panels show the diagonal Ci(k, k) of the full 3D power spectra 
for the full calculation (black lines) and for the Limber approximation (red, dark grey lines). The bottom panels show the modulus of 
the difference in the diagonal elements. This is for a Euclid wide survey (Refregier et al., 2010). 
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Figure 5. The Fisher matrix 2-parameter 1-cr constraints for the full 3D cosmic shear calculation and using the Limber approximation. 
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2010). 
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tends to too high a value of k, where theoretical uncertainties 
become a potential source of systematic error. None of this 
is a fundamental problem for tomography; it simply requires 
that the £ range chosen should be redshift-dependent - in- 
creasing £ max = r[z]fc max for the distant shells, and reducing 
it for nearby shells. 

In a similar manner to Section 2 we can also keep the sum- 
mation over galaxies explicit in the derivation of the tomo- 
graphic power spectrum. In this case we have a very similar 
expression to equation (12), with the kernel functions are 
now replaced by 



W iE 



i(r,n) 



h 



dr 



dr'p h (r'\r h ) 



^2ph(r'\r 



r' 



F k (r',r),(13) 



where we have a summed over the set of galaxies h — {g : 
n — Ar < r < ri + Ar} within a bin defined by a width 2Ar 
centred on n. Ph(r) is the redshift probability distribution 
for an individual galaxy within that bin. 



4.1 Convergence of 3D shear 



Here we investigate the predicted cosmological error con- 
straints from 3D cosmic shear as more modes are included. 
We assume a survey configuration as given in Section 2, 
and restrict the analysis to the regime £ max = 5000 and 
fcmax = 1.5hMpc _1 (with iV mo dos = 1000 k-modes evenly 
distributed over the range), we use the full posterior distri- 
butions from Bordoloi et al. (2009). 



the fiducial survey in this article is Ar survcy ~ 3000h _1 Mpc 
(where n ~ 3000 and r 2 ~ 0) such that Af mo d os ~ 
4500/1. For t = 150 this should result in convergence at 
■Wmodes, survey ~ 30. The photometric/n(z) smoothing scale 
typically occurs at Ar pn otoz ~ 150h _1 Mpc in comoving co- 
ordinates in this case we see that iV mo dos,photoz ~ 95, 000/£ 
from equation (14), for £ — 150 this should create conver- 
gence at Af mo dos ~ 600. In Figure 7 we see that for all 
£ ranges the convergence occurs at the photometric/n(z) 
smoothing scale. For high-i? modes the survey convergence 
is not visible since -/V mo d os , survey < 10- F° r l° w_ ^ modes the 
difference in the convergence regimes (30 and 600 k-bins) is 
clearly visible. 

Note that this simple illustration is approximate since the 
radial smoothing scale is also a function of the lensing ker- 
nel. In Figure 7 we also show the the FoM convergence for 
cr z (z)/(l + z) — 0.01 and 0.03. In both cases the the con- 
straints improve sharply over the range 1-20 modes. This 
suggests that below a z (z)/(l + z) ~ 0.03 the radial smooth- 
ing scale is not significantly reduced, due to the lensing ker- 
nel. The change in <J z {z) has a sub-dominant effect on the 
FoM especially at low numbers of modes, which is in agree- 
ment with tomographic studies (e.g. Bridle & King, 2007) 
and in the 3D limit with Heavens et al. (2006). 

We find that under the same survey assumptions, and a re- 
striction to modes £ < 5000, k < 1.5hMpc _1 that the 3D 
cosmic shear predictions agree very well with publicly avail- 
able tomographic code (for example iCosmo; Refregier et 
al., 2008) and other tomographic studies (e.g. Amara & Re- 
fregier, 2007; Bernstein, 2009). We note however that typical 
tomographic studies use highly non-linear modes, £ ~ 10 4 
and k ^> 1.5hMpc _1 , and consequently report tighter con- 
straints. 



Figure 6 shows how the Fisher matrix errors vary, for wo 
and )ii as well as the dark energy Figure of Merit (FoM= 

i/ V F ^on, F ^w a - (F^ Wa ) 2 ; Albrecht et al., 2006), as a 
function of the number of k modes sampled. We also show 
the 3D cosmic shear constraints using the Limber approxi- 
mation. In Figure 6 we also show how the constraints in the 
(wo,w a ) plane change as the number of modes is increased. 
We find that the cosmological constraints do not converge 
to the 3D (Limber) limit until the number of modes is > 800 
for low-£ modes. 



We emphasise here that our results are consistent and com- 
plimentary with tomographic studies. We find a convergence 
at the photometric smoothing scale. This is consistent with 
tomographic studies that find convergence at around 20 
redshift bins which is also consistent with the photometric 



smoothing scale where Ar « r s 



//^pho 



3000/150 = 



20. We recommend that in any analysis both tomography 
and 3D cosmic shear should be peformed 



We can understand the convergence if we consider that a 
separation of modes in fc-space Afc for a given £ corresponds 
to physical seperations Ar in the following way 



• Tomography allows one to explicitly scrutinise the red- 
shift dependence of the signal. 



JVmodcs nr 2 



(14) • 3D cosmic shear allows for an explicit radial and az- 

imuthal scale dependence to be investigated. 



where Ar = r± — r 2 and n — £/ki- We can use this as an ap- 
proximate model to investigate the convergence properties. 



We expect that convergence will occur in two regimes i) as 
radial modes enter the survey volume, ii) as radial modes 
become correlated due to the photometric/n(z) smoothing 
scale. To illustrate this we use equation (14). The depth of 



The redshift and scale dependence should be explicitly in- 
vestigated by both methods respectively. 
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Figure 6. Mode convergence of 3D cosmic shear. For the 3D calculation we use the Limber approximation, equation (11). The left panel 
shows how the inverse errors (l/c) for wo (middle green line), w a (lower red line) as well as the dark energy Figure of Merit (upper black 
line) (Albrecht et al., 2006) vary with the number of k modes, for for cr z (z)/(l + z) = 0.03 . The dotted horizontal lines show the 3D 
cosmic shear asymptotic predictions using the Limber approximation. The upper dashed (red) line shows the 3D cosmic shear prediction 
using the full (non-Limber) approximation. The inset shows the change in FoM over the range of 2 to 20 modes. The right panel shows 
the predicted 1-<t two-parameter constraints in the (wo,w a ) plane, the colour represents the number of k-mode bins that we define in 
the colour bar, the central white ellipse shows the Limber approximation 3D cosmic shear constraints. This is for a Euclid wide survey 
(Refregier et al., 2010). We use £ max = 5000 and fc max = l^hMpc" 1 . 



5 CONCLUSION 

In this article we have developed the 3D cosmic shear tech- 
nique in a number of ways. We have simplified the analysis 
using the Limber approximation, we have shown how to in- 
clude individual posterior redshift probability distributions 
from photo-zs in the analysis. The technique presented here 
removes the need for any intermediate p(z\z p ) estimation, 
allowing the estimator itself to become a direct function of 
the individual galaxy redshift probabilities. 

Finally, we have clarified the relationship between 3D cosmic 
shear and tomography, and demonstrated that tomography 
essentially provides an I- and shell-dependent sampling of 
physical wavenumber modes. 

To study individual redshift errors, we used a mock cata- 
logue of photometric redshift probabilities p(z) from Bor- 
doloi et al. (2009). We find that the cosmological errors can 
be mis-estimated by 10-50% by using an average p(z\z p ) 
in the estimator, but that biases on the cosmological pa- 
rameters are negligible. In order to calculate the biases we 
have generalised the bias formalism introduced in Knox et 
al. (1998) to the case where the cosmological signal is in the 
covariance not the mean. 

We have not investigated systematic effects that may be 
present in photometric redshift estimates. However we note 
that the outlier mitigation techniques and likelihood calibra- 
tion methodology of Bordoloi et al. (2009) is relevant and 
should be used in conjunction with the techniques presented 
here, and leave this for future investigation. In addition we 
refer to Kitching et al. (2008) (systematic effects on 3D cos- 
mic shear) where a factor of 2 degradation in FoM is ex- 
pected as a result of the primary weak lensing systematics. 



Using the extended Limber approximation of LoVerde & 
Afshordi (2008) we have found a much simpler expression 
for the 3D shear field in the high-^ regime. For angular 
wavenumbers £ > 100 the Limber approximation is essen- 
tially exact with a difference between the Limber and full 
calculation power spectra of < 10~ 4 ; a fractional difference 
of<10" 2 . 

Finally we explicitly derive the weak lensing tomographic 
power spectrum from the 3D cosmic shear field. We find 
that tomography essentially probes discrete sets of physical 
wavenumbers, which depend both on £ and on the distance 
to each individual tomographic bin. 

The choice of tomographic bin positions constrains the 3D 
modes which can be analysed. The modes which should be 
probed have to satisfy several constraints: first, the max- 
imum physical wavenumber should not be too high, to 
avoid the uncertain highly non-linear regime; secondly, as 
many modes as possible should be included, to reduce er- 
rors; thirdly the effective distance probed by adjacent modes 
should not be much less than that set by the photo-z error. 
The last of these is most easily effected with tomography, 
and the first two within 3D cosmic shear. 3D cosmic shear 
automatically deals (in principle) with the last constraint, 
through correlation of modes, but too finely-spaced modes 
may lead to near-singular covariance matrices. To deal with 
the first two constraints in tomography, the £ range chosen 
should be varied with shell. 

The advantages of 3D cosmic shear are in being able to 
ensure that appropriate physical modes are analysed, in an 
integrated analysis of the entire sample. On the other hand 
tomography allows for a straightforward investigation of the 
redshift-dependence of the shear signal. Consequently we 
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Figure 7. Left : We show the dark energy Figure of Merit (Albrecht et al., 2006) as a function of the number of k modes analysed. The 
lower/upper (red, green) lines are for a z (z)/(l + z) = 0.03 and 0.01 The inset shows the change in FoM over the range of 2 to 20 bins for 
<T z (z)/(l + z) = 0.03 (red, lower) and <r z (z)/(l + z) = 0.01 (green, upper). This is for a Euclid wide survey (Refregier et al., 2010). We use 
a { max = 5000 and A; max = 1.5hMpc . Right : For ff z (z) = 0.03 we show how the FoM changes for particular £ ranges (labelled near 
each line). Note that the in the full calculation, Figure 6, we use a range £ = 10—5000, and that the largest contribution to the Fisher 
matrix Fij = j AUFg^j (equation 3) is in the range i K, 1000 (Heavens et al., 2006). 



conclude that both approaches have their advantages and it 
is sensible to do both. 
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APPENDIX A: 3D COSMIC SHEAR 

In this appendix we generalize the results of Heavens et al., (2006) and Kitching et al. (2007). We compute the covariance of the 
harmonic shear coefficients utilizing the individual posterior galaxy redshift probability distributions, rather than simplified 
sample properties (i.e. an averaged p(z\z p )). 

Starting with equation (1) we can begin to derive a theoretical expression for the expected coefficient values by linking the 
shear values to the Newtonian potential $ via 7 = ^dd(f> where <j> is the lensing potential (d — di + idj) we can write a 
theoretical estimator like 

fy(M)=feyf#«)e- rf -°V(rg)j , (15) 

where i,j = {1,2} and 71 (M) = 3 (fn - f 22 ) and */ 2 {k,£) = f 12 = f 2 i. 
The covariance of equation (15) can be written as 



<f (*i,*)f (fe,/ a )> = J2Y1 Uj ^i(*i' , g)^(fe"-o)<^*>e- <€l - e +i£20 W{r 9 )W{r h ) (16) 

3 h 

where the covariance of the lensing potential can be written in terms of the Newtonian potential $ (note we assume General 
Relativity and no anisotropic stress, so the Newtonian and curvature potentials are identical): 

(<t>{rl6 9 W{rtO h )) = \ f dr' / dr"F K (r 3 ,r')F K (r h ,r")($(r\ 9 )<£>*(r", h )) (17) 

c Jo Jo 

where Fk(r,r') = Sk('i — r')/[Sk(r)Sk(r')] (Sk — sinh(fc), fc,sin(fc) for open, flat or closed geometries so Fk(r,r') = (r — r')/rr' 
for a flat universe). The real space Newtonian potential can also be written in terms of its spherical harmonic transform so 
that 

$( r , 9 )=(l) V2 J -^ J dkkMkr s Mk, I- ty LQ3 , (18) 

where t expresses the time-dependence. From now on we will identify this time dependence by an r dependence. Using Poisson's 
equation we can write the Newtonian potential's covariance (diagonal in t because of isotropy) in terms of the matter power 
spectrum P(k; r) as 

where A — 3ft m H 2 /2. Note that, as shown in Castro et al. (2003) we take advantage of an algebraic convenience by using the 
geometric mean of the matter power spectra, this is well justified since at large separations (where the approximation may 
break down) any pair correlations are down weighted by the form of the Bessel functions. 

This is fed into equations (15) then (16) so that we have an expression for the expected covariance of the 3D cosmic shear 
estimator 

2 



{T ij {k 1 ,e)f* ij {k2,i)) = amU* (J) (^) A 2 J2J2^ kir o^ i( - k2r ^ w ^ w ^) 



a h 

At' f " dr"F K (r 9 ,r')F K (r h ,r") [% \ j e (k'r')Mk'r")^/P(k';r')P(k';r"), (20) 

Jo J k a(r')a(r") 

the prefactor Af2, the angular size or area of the survey, comes from final integrations over angle (see Kitching et al., 2007 
for more information). 

We use this final expression in Section 2. 
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APPENDIX B: TOMOGRAPHY FROM 3D COSMIC SHEAR 

Here we derive the tomographic shear power spectra from the 3D cosmic shear covariance (equations 5, 6), assuming a global 

p(z\z p ). 

We start with the full 3D cosmic shear power spectrum as a function of azimuthal £ and two radial k wavenumbers. This is 
an expansion of the equations (5) and (6): 

Ci D {ki,k 2 ) = A 2 / dr g r 2 .n(r g )je(kir g ) / dr h rln{r h )je{k2r h ) / dr'p(r'\r g ) / dr"p(r"\r h ) 



(21) 



d ^ d "'' F %y' )F %'P I %Mk'f') M k'nP^;f')P^(k'; r") 



k' 2 ' 



where we see that four Bessel functions enter the expression. The prefactor A 2 is given by A 2 — £ 4 ■£ ^ 9Q 2 n H 4 /4. The 
p(r[z]) = p(r[z]|r[z p ]) is the redshift probability distribution, equivalent to p(z\z p ) in equation (6), and we condense the 
notation here to match the literature for the tomographic case. 

We now replace each Bessel function with its Limber approximated form (see equation 7) and find that we can rewrite the 
3D power spectrum 



^30, Limber/, ,x A 2 ( \ I 7 * 



^PWlll), 

a 2 (f) 



dff 



n(r 1 )r 1 n(r 2 )r2 I dr'p(r')F K (r',r) / dr"p(r")F K (r" , f) 



(22) 



where we have introduced two radial distances r\ and r 2 that are defined through the Limber approximation as r\ — l/ki 
and r 2 = l/k 2 . We can rewrite this in a more familiar form 



^i-iU, Limber / 7 7 \ 

C e (ki,k2) = 



9£l 2 m H 4 



dr 



P(£/r;r) W(n,r)W(r 2 ,r) 



4c 2 J a 2 (r) 

where we define a weight factor as 



W(ri,r) — [n(ri)r 2 ]r / dr'p(r'|n)^- 7 



, r — r 
r> 



(23) 



(24) 



We have written the geometric term in the flat geometry case for clarity FK=o{r,r') — (r — r')/rr' , and the r 2 term cancels; 
to match convention we add a factor of r to the weight factor that must be cancelled by a r 2 denominator. 

We note here that this expression is still a full 3D estimator in the Limber approximation. We use this expression in Section 
4 to derive the tomographic power spectra and relate this to the 3D cosmic shear power spectra. 



APPENDIX C: BIAS FOR SIGNAL IN THE COVARIANCE 

The bias formalism presented in Knox et al. (1998), Amara & Refregier (2007) and Taylor et al. (2007) assumes that the 
parameter dependency comes through the mean of the signal. Here we show how the bias formalism can be extended to the 
case where the signal is in the covariance. 

We start with a general expansion, in A0, of the (log) likelihood about the maximum likelihood position 6» + A8 



£(0o + AG) = £(6> ) + A0idiC(6 o ) + . 



(25) 



where Go in this case is the old position of the maximum likelihood and AG is the bias in parameters to a new maximum 
likelihood. By taking the derivative with respect to the parameters we have 



djC(G + AG) = 2djC{G ) + A6Ad 3 C{6 {) ) + . . . = 



(26) 



where the log-likelihood is maximised at £(Go + AG). By taking the expectation value of this and rearranging we find that 
the bias bi in the i th parameter can be written as 



h^A6 l ^Fr j 1 (2d J C{G Q )), 



(27) 
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in a completely general way. Note that this makes the assumption that the Fisher matrix is unaffected to first order (this is 
a good approximation for small systematics as shown by Joachimi & Bridle, 2009). 

By taking the result from Tegmark, Taylor & Heavens (f997) in the Gaussian case we can write the first derivative of the 
log-likelihood as 

(2£,i ) t = Tr [Ci x C Ul -Ci x C Ul C^D + C e D„ ] (28) 

as a function of £, where we replace di —,i. D is the data matrix {x — /j,)(x — /*)'. Here and for the remainder of this Appendix 
we assume an all sky survey such that the covariance Ci is diagonal in £ (see Heavens et al., 2006; Kitching et al., 2007). 
We now assume that the mean is not a function of parameters and only consider the case where the parameters are in the 
covariance D,i = 0; in the case that the signal is in the mean and Ce,i = we recover the result of Amara & Refregier (2007) 
and Joachimi & Bridle (2009). Note that the quantity in the square brackets is a {£, k, k') 3D matrix, and that the derivative 
of the log-likelihood is a function of t. 

To include a systematic we rewrite the covariance as some true covariance plus a systematic C — C l + C s . By taking the 
expectation value of equation (28) we have 

{2C ti ) t = Tr[(C| + Ciy x CU -{C\ + Ciy'Cli {Cl + C|)- 1 (D)](2£, l ) t = Tr [{C\ + Ct^CU (i - {C\ + Ciy'd)] (29) 

where the expectation of the data is the true covariance (D) = C\. Finally we note that the total log- likelihood is in general 
a sum over £ — (t x ,ty) modes C — y^ . , £e x ,e y , so a similar summation must be performed over the derivatives of the 
log- likelihood. For the Fisher matrix and for the bias here we write this summation as an integral over the modulus t — \t\, 
and account for the density of state in (£ x ,ty) (see Appendix B in Kitching, Taylor, Heavens, 2007). 

We can now write an expression for the bias caused by a systematic function in the case of the parameters being in the 
covariance like 

bi = Fr/{g { A4>, I d^Tr [{CJ + C?)" 1 ^- (/ - {C\ + C?)" 1 ^)] }. (30) 

In this final step we have also added an integration over £-space where g is the density of states in £; this is exactly the same 
integration that is performed when calculating the Fisher matrix (see Appendix B in Kitching, Taylor, Heavens, 2007) . If the 
systematic is very low \C S \ <g; C then (2£,; ) = and the bias is effectively zero. 
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